Microbial individuality: how single-cell heterogeneity enables population level strategies.

Much of our knowledge of microbial life is only a description of average population behaviours, but modern technologies provide a more inclusive view and reveal that microbes also have individuality. It is now acknowledged that isogenic cell-to-cell heterogeneity is common across organisms and across different biological processes. This heterogeneity can be regulated and functional, rather than just reflecting tolerance to noisy biochemistry. Here, we review recent advances in our understanding of microbial heterogeneity, with an emphasis on the pervasiveness of heterogeneity, the mechanisms that sustain it, and how heterogeneity enables collective function.The research has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement 338060. The work in the Locke laboratory is also supported by a fellowship from the Gatsby Foundation (GAT3273/GLC).This is the author accepted manuscript. The final version is available from Elsevier via http://dx.doi.org/10.1016/j.mib.2015.01.00

Widespread inter-individual gene expression variability in Arabidopsis thaliana.

A fundamental question in biology is how gene expression is regulated to give rise to a phenotype. However, transcriptional variability is rarely considered although it could influence the relationship between genotype and phenotype. It is known in unicellular organisms that gene expression is often noisy rather than uniform, and this has been proposed to be beneficial when environmental conditions are unpredictable. However, little is known about inter-individual transcriptional variability in multicellular organisms. Using transcriptomic approaches, we analysed gene expression variability between individual Arabidopsis thaliana plants growing in identical conditions over a 24-h time course. We identified hundreds of genes that exhibit high inter-individual variability and found that many are involved in environmental responses, with different classes of genes variable between the day and night. We also identified factors that might facilitate gene expression variability, such as gene length, the number of transcription factors regulating the genes and the chromatin environment. These results shed new light on the impact of transcriptional variability in gene expression regulation in plants

Co-expression Networks From Gene Expression Variability Between Genetically Identical Seedlings Can Reveal Novel Regulatory Relationships.

Co-expression networks are a powerful tool to understand gene regulation. They have been used to identify new regulation and function of genes involved in plant development and their response to the environment. Up to now, co-expression networks have been inferred using transcriptomes generated on plants experiencing genetic or environmental perturbation, or from expression time series. We propose a new approach by showing that co-expression networks can be constructed in the absence of genetic and environmental perturbation, for plants at the same developmental stage. For this, we used transcriptomes that were generated from genetically identical individual plants that were grown under the same conditions and for the same amount of time. Twelve time points were used to cover the 24-h light/dark cycle. We used variability in gene expression between individual plants of the same time point to infer a co-expression network. We show that this network is biologically relevant and use it to suggest new gene functions and to identify new targets for the transcriptional regulators GI, PIF4, and PRR5. Moreover, we find different co-regulation in this network based on changes in expression between individual plants, compared to the usual approach requiring environmental perturbation. Our work shows that gene co-expression networks can be identified using variability in gene expression between individual plants, without the need for genetic or environmental perturbations. It will allow further exploration of gene regulation in contexts with subtle differences between plants, which could be closer to what individual plants in a population might face in the wild

Frequency doubling in the cyanobacterial circadian clock

Organisms use circadian clocks to generate 24-h rhythms in gene expression. However, the clock can interact with other pathways to generate shorter period oscillations. It remains unclear how these different frequencies are generated. Here, we examine this problem by studying the coupling of the clock to the alternative sigma factor $\textit{sigC}$ in the cyanobacterium $\textit{Synechococcus elongatus}$. Using single-cell microscopy, we find that $\textit{psbAI}$, a key photosynthesis gene regulated by both $\textit{sigC}$ and the clock, is activated with two peaks of gene expression every circadian cycle under constant low light. This two-peak oscillation is dependent on $\textit{sigC}$, without which $\textit{psbAI}$ rhythms revert to one oscillatory peak per day. We also observe two circadian peaks of elongation rate, which are dependent on $\textit{sigC}$, suggesting a role for the frequency doubling in modulating growth. We propose that the two-peak rhythm in $\textit{psbAI}$ expression is generated by an incoherent feedforward loop between the clock, $\textit{sigC}$ and $\textit{psbAI}$. Modelling and experiments suggest that this could be a general network motif to allow frequency doubling of outputs.This research was made possible by the award of a European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement 338060. The work in the Locke laboratory is further supported by a fellowship from the Gatsby Foundation (GAT3272/GLC) and a Fellowship from the Human Frontier Science Program (CDA00068/2012)

Coordinated circadian timing through the integration of local inputs in Arabidopsis thaliana.

Individual plant cells have a genetic circuit, the circadian clock, that times key processes to the day-night cycle. These clocks are aligned to the day-night cycle by multiple environmental signals that vary across the plant. How does the plant integrate clock rhythms, both within and between organs, to ensure coordinated timing? To address this question, we examined the clock at the sub-tissue level across Arabidopsis thaliana seedlings under multiple environmental conditions and genetic backgrounds. Our results show that the clock runs at different speeds (periods) in each organ, which causes the clock to peak at different times across the plant in both constant environmental conditions and light-dark (LD) cycles. Closer examination reveals that spatial waves of clock gene expression propagate both within and between organs. Using a combination of modeling and experiment, we reveal that these spatial waves are the result of the period differences between organs and local coupling, rather than long-distance signaling. With further experiments we show that the endogenous period differences, and thus the spatial waves, can be generated by the organ specificity of inputs into the clock. We demonstrate this by modulating periods using light and metabolic signals, as well as with genetic perturbations. Our results reveal that plant clocks can be set locally by organ-specific inputs but coordinated globally via spatial waves of clock gene expression

Interpretation of morphogen gradients by a synthetic bistable circuit.

During development, cells gain positional information through the interpretation of dynamic morphogen gradients. A proposed mechanism for interpreting opposing morphogen gradients is mutual inhibition of downstream transcription factors, but isolating the role of this specific motif within a natural network remains a challenge. Here, we engineer a synthetic morphogen-induced mutual inhibition circuit in E. coli populations and show that mutual inhibition alone is sufficient to produce stable domains of gene expression in response to dynamic morphogen gradients, provided the spatial average of the morphogens falls within the region of bistability at the single cell level. When we add sender devices, the resulting patterning circuit produces theoretically predicted self-organised gene expression domains in response to a single gradient. We develop computational models of our synthetic circuits parameterised to timecourse fluorescence data, providing both a theoretical and experimental framework for engineering morphogen-induced spatial patterning in cell populations

Escherichia coli can survive stress by noisy growth modulation.

Gene expression can be noisy, as can the growth of single cells. Such cell-to-cell variation has been implicated in survival strategies for bacterial populations. However, it remains unclear how single cells couple gene expression with growth to implement these strategies. Here, we show how noisy expression of a key stress-response regulator, RpoS, allows E. coli to modulate its growth dynamics to survive future adverse environments. We reveal a dynamic positive feedback loop between RpoS and growth rate that produces multi-generation RpoS pulses. We do so experimentally using single-cell, time-lapse microscopy and microfluidics and theoretically with a stochastic model. Next, we demonstrate that E. coli prepares for sudden stress by entering prolonged periods of slow growth mediated by RpoS. This dynamic phenotype is captured by the RpoS-growth feedback model. Our synthesis of noisy gene expression, growth, and survival paves the way for further exploration of functional phenotypic variability

Global parameter search reveals design principles of the mammalian circadian clock

Background: Virtually all living organisms have evolved a circadian (~24 hour) clock that controls physiological and behavioural processes with exquisite precision throughout the day/night cycle. The suprachiasmatic nucleus (SCN), which generates these ~24 h rhythms in mammals, consists of several thousand neurons. Each neuron contains a gene-regulatory network generating molecular oscillations, and the individual neuron oscillations are synchronised by intercellular coupling, presumably via neurotransmitters. Although this basic mechanism is currently accepted and has been recapitulated in mathematical models, several fundamental questions about the design principles of the SCN remain little understood. For example, a remarkable property of the SCN is that the phase of the SCN rhythm resets rapidly after a 'jet lag' type experiment, i.e. when the light/ dark (LD) cycle is abruptly advanced or delayed by several hours. Results: Here, we describe an extensive parameter optimization of a previously constructed simplified model of the SCN in order to further understand its design principles. By examining the top 50 solutions from the parameter optimization, we show that the neurotransmitters' role in generating the molecular circadian rhythms is extremely important. In addition, we show that when a neurotransmitter drives the rhythm of a system of coupled damped oscillators, it exhibits very robust synchronization and is much more easily entrained to light/dark cycles. We were also able to recreate in our simulations the fast rhythm resetting seen after a 'jet lag' type experiment. Conclusion: Our work shows that a careful exploration of parameter space for even an extremely simplified model of the mammalian clock can reveal unexpected behaviours and non-trivial predictions. Our results suggest that the neurotransmitter feedback loop plays a crucial role in the robustness and phase resetting properties of the mammalian clock, even at the single neuron level

Transancestral fine-mapping of four type 2 diabetes susceptibility loci highlights potential causal regulatory mechanisms.

To gain insight into potential regulatory mechanisms through which the effects of variants at four established type 2 diabetes (T2D) susceptibility loci (CDKAL1, CDKN2A-B, IGF2BP2 and KCNQ1) are mediated, we undertook transancestral fine-mapping in 22 086 cases and 42 539 controls of East Asian, European, South Asian, African American and Mexican American descent. Through high-density imputation and conditional analyses, we identified seven distinct association signals at these four loci, each with allelic effects on T2D susceptibility that were homogenous across ancestry groups. By leveraging differences in the structure of linkage disequilibrium between diverse populations, and increased sample size, we localised the variants most likely to drive each distinct association signal. We demonstrated that integration of these genetic fine-mapping data with genomic annotation can highlight potential causal regulatory elements in T2D-relevant tissues. These analyses provide insight into the mechanisms through which T2D association signals are mediated, and suggest future routes to understanding the biology of specific disease susceptibility loci

Genetic fine mapping and genomic annotation defines causal mechanisms at type 2 diabetes susceptibility loci.

We performed fine mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in or near KCNQ1. 'Credible sets' of the variants most likely to drive each distinct signal mapped predominantly to noncoding sequence, implying that association with T2D is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine mapping implicated rs10830963 as driving T2D association. We confirmed that the T2D risk allele for this SNP increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease